Schwinger-Keldysh semionic approach for quantum spin systems 
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We derive a path-integral Schwinger-Keldysh approach for quantum spin systems. This is achieved 
by means of a semionic representation of spins as fermions with imaginary chemical potential. The 
major simplifying feature in comparison with other representations (Holstein-Primakoff, Dyson- 
Maleev, slave bosons/fermions etc) is that the local constraint is taken into account exactly. As a 
result, the standard diagram technique with usual Feynman codex is constructed. We illustrate the 
application of this technique for the Neel and spin-liquid states of the AFM Heisenberg model. 
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Since a long time physicists are aware of the fact 
that spin operators which commute on different sites 
and anticommute on the same site are neither Fermi 
nor Bose operators. Less convergent opinions exist on 
whether fermionizations or bosonizations or none of those 
should be used to take care of spin statistics in many 
body quantum theory. At least the answers appear to 
be linked to the kind of physical problem considered. 
Widely accepted is the view that path integral represen- 
tations and diagrammatic expansions for spin systems 
are thus substantially more complicated than those of 
pure fermion/boson systems. Many variants of diagram 
technique j^] that are based on different representation 
of spins such as Bose 0]- [Q, Fermi Js]- [^|, Majorana 
supersymmetric ||] or Hubbard [|^| operators have 
been proposed. Another method to treat spin Hamilto- 
nians is based on direct representation of coherent states 
for spins (nonlinear cr model, see e.g ^). Some of these 
techniques being applicable only at low temper- 

atures or in large spin (S'3>1) limit, describe nevertheless 
well excitations in ordered magnets (ferromagnetic and 
antiferromagnetic magnons), but fail to provide rigorous 
calculations in strongly correlated systems such as Kondo 
lattices or quantum magnets. Other techniques, based on 
a successful choice for the hierarchy of coupling constants 
are mainly restricted to equilibrium situations. The fun- 
damental problem which is at the heart of the difficulty 
is the local constraint. On one hand, any representa- 
tion of spin operators as bilinear combination of Fermi 
or Bose quasiparticles makes the dimensionality of the 
Hilbert space, where these operators act, greater than 
the dimensionality of the Hilbert space for spin opera- 
tors. As a result, the spurious unphysical states should 
be excluded from the consideration resulting in substan- 
tial complication of corresponding rules of diagrammati- 
cal summation. On the other hand, there is no Wick the- 
orem directly for spin operators but the Gaudin theorem 
instead (see also |7[|lO|). It can however not avoid 
complications in diagram techniques based on Hubbard 
operators, rendering the resummation of diagram series 
in many cases practically uncontrollable. The exclusion 
of double-occupied and empty states for impurity spins 



interacting with a conduction electron bath (single im- 
purity Kondo problem) ||^ is cured by an infinite chemi- 
cal potential for Abrikosov pseudofermions. It works for 
dilute spin subsystems, where all spins can be consid- 
ered independently. Unfortunately, attempts to gener- 
alize this technique to the lattice of spins result in the 
replacement of the local constraint (the number of parti- 
cles on each site is fixed) by a so-called global constraint 
(in the saddle point approximation), where the number 
of particles is fixed only as an average value for the whole 
crystal. There is no reason to believe that such an ap- 
proximation is a good starting point for the description 
of strongly correlated systems. Besides, it is very diffi- 
cult to take into account the fluctuations related to the 
replacement of a local constraint by a global one. 

An alternative approach for spin Hamiltonians free of 
the local constraint problem has been proposed in the 
pioneering paper of Popov and Fedotov (PF) |lj]. Based 
on exact fermionic representation for 5=1/2 and 5=1 
operators, where fermions are treated as quasiparticles 
with imaginary chemical potential, these authors demon- 
strated the power and simplification of the correspond- 
ing Matsubara diagram technique. For these two special 
cases the Matsubara frequencies are cj„=27rT(n-|-l/4) for 
5=1/2 and uJn—2nT{n + 1/3) for 5=1 providing a rig- 
orous description of (and restricted to) the equilibrium 
situation. The semionic representation used by PF is nei- 
ther fermionic, nor bosonic, but reflects the fundamental 
Pauli nature of spins. Later the generalization of the PF 
technique for arbitrary spin [ p^ was derived by intro- 
ducing proper chemical potentials for spin fermions. The 
goal of this paper is to derive a method for nonequilib- 
rium systems, which allows to treat quantum spin Hamil- 
tonians on the same footing as Fermi or Bose systems. 

A long time ago Keldysh jlj] proposed a novel ap- 
proach for the description of kinetic phenomena in met- 
als. This approach was found especially fruitful for nor- 
mal metals and, in many recent applications, for su- 
perconductors ||l6|| , for disordered interacting (normal or 
superconducting) electron liquids for example. The 
previous application of the real-time formalism to the 
quantum theory of Bose-Einstein condensation (BEG) 
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p8[ allowed the derivation of a Fokker-Planck equation, 
which describes both kinetic and coherent stages of BEC. 
Moreover developed the closed-time path integral 
formalism for aging effects in quantum disordered sys- 
tems being in contact with an environment. The Keldysh 
formalism in application to disordered systems (see |2^] 
- [^) also attracted interest some time ago as an al- 
ternative approach to the replica technique. The main 
advantage of closed-time contour calculations is an auto- 
matic normalization (disorder independent) of the parti- 
tion function. In this paper we derive the Keldysh for- 
malism for quantum spin systems (e.g. Heisenberg clean 
and disordered magnets, Kondo lattices), which is based 
on Popov-Fedotov ideas of semionic representation. 

We reformulate the PF-concept adopting it to real- 
time formalism. As an example we consider S = 1/2. 
As it was first shown in [|l2|, the partition function of 
a spin system with Hamiltonian Hs can be replaced by 
the partition function of an effective "fermionic" system 
with Hamiltonian Hp as follows 



(1) 



where (3=1 /T and usual "fermionic" representation of 
spin similar to e.g. Abrikosov pseudofermions |^] 
is used: 5+=/|/x, S-=flf^, s^ = \{f\f^-flf^), and 

Representing spins as bilinear combinations of Fermi 
operators, we enlarged by a factor of two the Hilbert 
space of the Hamiltonian. In addition to physical 
states |1,0) and |0, 1) two unphysical states |1,1) and 
|0, 0) are introduced. Nevertheless in the average over 
all states unphysical states cancel each other, since 

TrunphysieMTiT^I'^))^^ = (T«)° + {Ti? = 0. This rep- 
resentation being of semionic origin results in conven- 
tional Matsubara diagram technique with ujn = 2TTT{n + 
1/4) or u)n = 27rr(n -I- 3/4) depending on the sign in 
expressions (|]). Besides, one can introduce the auxiliary 
distribution function for quasiparticles [E3[ 
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Jn-€ e±*'^/2exp(/3e) + 1 



(2) 



where signs ± in the exponent (|^) are the same as in 
(|l|). We note that, since auxiliary Fermi fields do not 
represent the true quasiparticles of the problem helping 
only to treat properly the spin operators, the distribu- 
tion function for these objects in general should not be 
a real function, e.g. /'■^/^•' = n(2e) =F |sech(e/T) where 
n{x) = [exp(x/r) -f 1]"^ is the standard Fermi distribu- 
tion function. As we shall see for arbitrary value of spin, 
1 — 2Ref^^\e) = Bs{t/T) is expressed in terms of the 
Brillouin function Bs{x) = (1 -f ^)coth((l -t- ■:^)x) — 
^coth(^), e.g. for S = 1/2, Bi/2{x) = tanh(a;). We 
also note that in the T — > 0-limit the imaginary part of 
satisfies the identity Imf^^^^^{x) = Ti'^TS{x)/2. 



The spin correlation functions of any order can be ex- 
pressed in terms of the two-component field i/'"^=(/t /i)' 

(5-(ii)...5-(t„)> = rr(po(^vr»,,...(V'Ve:v)tj 

where po = exp(— /3i/o)/rr cxp(— /3i7o) is the density 
matrix and a denotes Pauli matrices. We included the 
term inNp / (2/3) into the Hamiltonian Hq = —h ± 

ittT /2 Np'^ of noninteracting spins in a uniform ex- 
ternal magnetic field ft,, since it exists both in numerator 
and denominator of po- 
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FIG. 1. Keldysh double side contour going along real time 
axis — oo +00 — oo and "closed" in imaginary time. 

Following the standard route |Q we can express the 
partition function of the problem containing spin opera- 
tors as a path integral over Grassmann variables Tp^if) 



(3) 



where actions A and Aq are taken as an integral along 
the closed-time contour Ct + Cr which is shown in Fig.l. 
The contour is closed at i = ~oo + iT since exp(— /?i/o) = 
Tr exp ^— Hadrj . We denote the V fields on upper 

and lower sides of the contour Ct as tpi and 7/12 respec- 
tively. The fields ^E" stand for the contour Cr- These 
fields provide matching conditions for 'ipi^2 and are ex- 
cluded from final expressions. Taking into account the 
semionic boundary conditions for generalized Grassmann 
fields *'^(/3)=i*''(0), ^^'{f3) = -i^^'{0) one gets the 
matching conditions for 7/11.2 at i = ±00. 

7/;^(-oo) = zV'^(-oo)e'3'*<', Vi(+oo) = V2(+oo) (4) 

The correlation functions can be represented as func- 
tional derivation of the generating functional 



Z['ij]/Zo = J DtPDiP exp l^iA + i j) dt{fja^iP + ^Pa""!]) 

where i] represents sources and cr^ matrix stands for 
"causal" and "anti-causal" orderings along the contour. 

The on-site Green's functions (GF) which are matri- 
ces 4x4 with respect to both Keldysh (lower) and spin 
(upper) indices are given by 



After a standard shift-transformation [g4| of fields ijj the 
Keldysh GF of free PF-fermions assumes the form 
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G^(e) = G, 



l-fe -fe 



Gr. 
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where the retarded and advanced GF's are 

G^«'^)"(.) = (6 + a^V2±*^)-\ fe^f^'/'He). (5) 

The interdependence of matrix elements of the GF in 
Keldysh-space is more transparent after rotation 
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(6) 



where G^= -i27r5(e ± h/2) [B^/^ie/T) ± i sech(e/r)] . 
We emphasize, that unUke diagrammatic techniques 
for Fermi and Bose operators, the off-diagonal element 
(Keldysh-component) in semionic representation is ex- 
pressed in terms of Brillouin function, containing correct 
information about occupied states. We recall that di- 
agonal elements of the matrix (|^) in "triangular" repre- 
sentation satisfy the Dyson equation providing the exact 
description of the system. The equation of motion for 
G^ generally constitutes the quantum-kinetic equation. 

Let us illustrate the application of the Schwinger- 
Keldysh formalism for spin Hamiltonians. We consider 
the Heisenberg model with nearest neighbour interaction 
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We firstly discuss the Neel solution for the Heisenberg 
model with isotropic antiferromagnetic (AFM) exchange 
(J < 0). Applying the PF-transformation to the parti- 
tion function one obtains the action as an integral along 
the closed-time Keldysh-contour 

A = Ao + A^nt = Ao + £dtY, J(q)5q(i)^-q(0 (7) 



where Aq corresponds to noninteracting fermions 



Aa 



We denote Jq=^X]<i>6*''', t^q=Jq/Jo and apply 8- 
component PF representation with ?A"^=('0k '/'k+q)' 
where Q = (tt, ...,7r) for hypercubic lattice. To decouple 
the four-fermion term along the Keldysh-contour with 
the help of Hubbard- Stratonovich transformation we in- 
troduce the two-Keldysh-component vector (Bose) field 
= ($1 $2)- As a result we obtain 



A^nt = --Tr($^ Jq V^$q) + Tr(^$a7^^) (9) 



Now we integrate out ip fields and express the effective 
action in terms of $ fields 

■^eff = -^rr($^ Jq V^$q) + Tr In (Gq ^ + $^a7^) 



where 7^=((t^±1)/2 acts in Keldysh space. Since in gen- 
eral $ is a time- and space-dependent fluctuating field 
the partition function (H) cannot be evaluated exactly. 
Nevertheless, when a magnetic instability occurs, we can 
represent the longitudinal component of this field as a su- 
perposition of a staggered time- independent part (" stag- 
gered condensate" ) and a fluctuating field 



$^(q, UJ) = AAjq7^,5q,Q^(cj) + (/.^(q, ij), 



(10) 



where A/" is a staggered magnetization and <i>^ (q, uj) — 
(f>^{q,uj) with the matching conditions at t = ±00 in- 
volving the dispersions of excitations ujp 

0±(-oo) = 0±(-oo)e^^-p, (/)f(+oo)-0±(+oo). (11) 

We expand Trln(GQ"^ + 4>^a^^) in accordance with 

Tr ln(...) - TrlnGo 1 + V (Go0^a7^)" (12) 




FIG. 2. Feynman diagrams contributing to dispersion (a) 
and damping (b) of magnons. Solid line denotes PF fermions. 

The spectrum of the excitations (AFM magnons) 
can be defined as poles of the transverse GF 
D+-=i:)(x,i)=-z(rc<?!)^(x,i)0i'(O,O)). The procedure 
of the calculation of this GF is similar to that for a 
"fermionic" GF. Introducing the sources and evaluating 
(O) one gets 



where the retarded and advanced magnons GF's are 

i?^'^(p,w) = {u;-Luip)±i6y\ = (exp(/3c^)-l)-i 

The magnon spectrum is determined by the zeros of the 
expression — HJ [uu) (see Fig. 2a) in equilibrium 



LOp = \JoW^T^l => c|p|, N = tanh (^^) (13) 

The magnon-damping is defined by four-magnon pro- 
cesses HJ ' , shown in Fig2.b. The derivation of the 
kinetic equation and calculation of magnon damping is 
reserved here for a detailed publication. The results ( p^ ) 
(and similar for quantum FM) are in full agreement with 
the spin- wave theory (see, for example, j|] and ^Q)- 

The second possibility to decouple the four-fermion 
term in the Heisenberg model is provided by the bi-local 
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scalar bosonic field A^- depending on two sites. Intro- 
ducing new coordinates R={Ri+Rj)/2, p=Ri—Rj and 
applying a Fourier transformation we obtain the effec- 
tive action 

A// - -\Tr{hl^^J^^_^yKp,,)+Tr\n{G^' - A^7^) 

This effective action describes the nonequilibrium quan- 
tum spin-liquid (SL). We confine ourselves to consider 
the uniform phase ||2^ of Resonant Valence Bonds. It is 
sutable to rewrite the functional in new variables, namely 
amplitude A and phase Q=pA{R)^ according to formula 

A<^^>(i?, p) = Aip}jrc^P {^pMR)) (14) 

The exponent in ( p^ stands for gauge fiuctuations to 
be taken in eikonal approximation. The spectrum of 
excitations in the uniform SL is defined by the zeros 
of 7r«^"'3=Tr(pV(G«+,Gf +Gf+,G^)+5„^/(JpA)) in 
equilibrium [g6| and is purely diffusive (see e.g. ||25|] ) 

t^ = zJA|q|3, A = -^Kq)tanh(':^') (15) 

The quantum kinetic equation for nonequilibrium spin 
liquids can be obtained by taking into account the higher 
order diagrams similarly to Fig. 2b with current-like ver- 
tices and will be presented elsewhere. 

We discuss finally the Schwinger-Keldysh formalism for 
spins S > 1/2. As shown by Popov and Fedotov for S'=l, 
it is possible to eliminate the unphysical states by intro- 
ducing three-component fermions i/'"^=(/t /o fi) with 
imaginary chemical potential A=— zttT/S. The boundary 
conditions for on the imaginary part of the contour Cr 
read as follows 4'^(/3)=e"/34'^(0), *'^(/3)=e-*'^/3^'^(0). 
As a result, the distribution function in equilibrium is 
/(i)(e)=l/[e±W3exp(e/T) + l] . Thus the Schwinger- 
Keldysh formalism with 6x6 matrices for GF (^) and 
fe=f^^\^) in equilibrium is obtained. The off-diagonal 
Keldysh component is given by 

- -z2^(5(e±/i)[i?i(e/r)±*V3sinh(^)/sinh(^)]. 

For arbitrary spin values iS'>l there is no unique 
imaginary chemical potential for 25* -1-1 component PF- 
fermions, but instead they are distributed on each lattice 
site j according to 

[^^-1/2] 2i f 2l + l\ 

where \i^iT:T{2l + l)/{2S + l) |3|. Thus, the Schwinger- 
Keldysh approach can be generalized for arbitrary spin 
values in the same fashion as for 5=1/2 and 5=1. 

Summarizing, we derived the technique applicable for 
nonequilibrium dynamics of quantum spin systems. Un- 
like other techniques this approach takes into account the 



constraint rigorously and allows one to treat spins on the 
same footing as Fermi and Bose systems. The method 
derived can be applied especially to problems where the 
local constraint becomes important, e.g. quantum phase 
transition in clean and disordered magnets, spin glasses, 
Kondo lattices, nonequilibrium Kondo systems etc. 
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